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We present a model for electron- and neutrino-scattering oflt nucleons and nuclei focussing on 
the quasielastic and resonance region. The lepton-nucleon reaction is described within a relativistic 
formalism that includes, besides quasielastic scattering, the excitation of 13 A''* and A resonances 
and a non-resonant single-pion background. Recent electron-scattering data is used for the state-of- 
the-art parametrizations of the vector form factors; the axial couplings are determined via PCAC 
and, in the case of the A resonance, the axial form factor is refitted using neutrino-scattering 
data. Scattering off nuclei is treated within the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) 
framework that takes into account various nuclear effects: the local density approximation for the 
nuclear ground state; mean-field potentials and in-medium spectral functions. Results for inclusive 
scattering off Oxygen are presented and, in the case of electron-induced reactions, compared to 
experimental data and other models. 

PACS numbers: 25.30.-c 



I. INTRODUCTION 

Scattering electrons off the nucleon, a wealth of infor- 
mation on its properties such as, e.g., structure functions, 
charge radius, the distribution of its quark and gluon con- 
stituents and its excitation spectrum has been gathered. 
Under investigation are also how these properties change 
when embedding the nucleon or the resonances in a nu- 
clear medium. Both the resonances and the nucleon ac- 
quire complex self-energies within in the medium which 
may lead both to mass-shifts and modifications of the life- 
times. Calculating these with a given interaction model 
and comparing such predictions to the measured quan- 
tities, one is testing our understanding of the hadronic 
many body problem. Such studies have first been per- 
formed via inclusive experiments with nuclear targets 
where only the outgoing electron was detected (for a re- 
cent review cf., e.g., Ref. fl). Studying semi-inclusive 
processes in the intermediate energy regime with an en- 
ergy transfer of 0.1 — 2 GeV, one gains further sensitivity 
on in-medium changes of baryonic resonances and the 
nucleon since such modifications may also lead to unex- 
pected final state interaction patterns. A related topic 
of present interest is the modification of mesons such as, 
e.g., a and ui Q within the medium due to chiral 
symmetry restoration and/or coUisional broadening. A 
correct understanding of final state interactions is neces- 
sary to distinguish profane (e.g., pion rescattering in the 
nucleus) from spectacular (e.g., chiral symmetry restora- 
tion) effects. 

Closely related is the scattering of neutrinos on nu- 
cleons, where current and future experiments propose to 
shed light on the nucleon axial form factor and its strange 
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quark content. However, all those neutrino experiments 
use nuclear targets and, therefore, all measured cross sec- 
tions incorporate nuclear effects [J] . To draw conclusions 
on the underlying vN process it is therefore necessary 
to understand these corrections. Nowadays, the inter- 
est in neutrino nucleus reactions is driven by the dis- 
covery of neutrino oscillations — there is an extensive ex- 
perimental effort aiming at a precise determination of 
neutrino oscillation parameters as mixing angles, neu- 
trino mass squared differences and possible CP violation 
[E Hj H B S ■ This demands for an equally precise 
knowledge of the neutrino nucleus interaction process. A 
critical quantity is the neutrino energy which can not 
be measured directly but has to be reconstructed from 
observables, i.e., in the case of charged current (CC) re- 
actions from the outgoing charged lepton, or in the case 
of neutral currents (NC) from the hadronic debris leav- 
ing the nucleus. Neutrino induced pion production is 
strongly influenced by nuclear effects, its understanding 
is crucial since NC tt^ production is a major background 
in h'e appearance experiments, while CC tt"*" production 
introduces a background to disappearance searches. A 
good knowledge of neutrino-nucleus interactions is thus 
necessary to minimize the systematic uncertainties in 
neutrino fluxes and backgrounds. 

Electron and neutrino scattering off the nucleon are 
closely interconnected and one can treat both processes 
within the same formalism. Furthermore, the proper de- 
scription of the electron nucleus interaction serves as a 
necessary benchmark for the nuclear corrections in the 
neutrino nucleus reaction. 

Many authors investigated electron or neutrino scat- 
tering on nucleons, but only a few study both within 
the same model. In the intermediate energy regime, the 
lepton-nucleon cross section is dominated by elastic scat- 
tering and pion production via resonance excitations or 
non-resonant processes. One approach to describe pion 
production is based on helicity amplitudes derived in a 
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relativistic quark model [ll| . Those have been adopted to 
describe both, electro- and neutrino production (via CC 
and NC) of resonances [13, ■ Nieves and collaborators 
extended their model for eN e'N'ir [l^l to describe 
CC and NC pion production on the nucleon [l^ . Besides 
the dominant A contribution they included background 
terms required by chiral symmetry. Also the model of 
Sato and Lee for neutrino induced pion production on the 
nucleon inclu ding A excitation and non-resonant back- 
ground terms [iq is based on their approach for electrons 

The influence of nuclear effects on both, electron and 
neutrino scattering cross sections, has been investigated 
by Benhar et al. using an impulse approximation model 
with realistic spectral functions obtained from nuclear 
many body theory calculations [H, [l^ . In particular in 
the quasielastic peak region, they achieve good agreement 
to inclusive electron scattering data. Nuclear effects in 
the QE region have also been investigated in detail by 
Nieves et al. for electrons 14 1 and neutrinos [13, [21| 



where they have included, among other nuclear correc- 
tions, long range nuclear correlations. Also this approach 
describes inclusive electron scattering data with impres- 
sive agreement. A relativistic Green's function approach 
was applied by Meucci et al. [13, to inclusive elec- 
tron as well as to inclusive neutrino nucleus reactions. 
Butkevich et al. [H, addresses both neutrino and 
electron scattering with special emphasis on the impact 
of different impulse approximation (lA) schemes: plane 
wave lA (PWIA) and relativistic distorted wave lA (RD- 
WIA). Also the Ghent group apply RPWIA and RDWIA 
models to neutrino and electron scattering in the QE re- 
gion 26]; l ately they extended their framework to pion 
production [27| . 

In this article, we concentrate on two issues: First, we 
consider electron and neutrino — both charged (CC) and 
neutral current (NC) — scattering off the nucleon. Then, 
we investigate the influence of the nuclear medium on in- 
clusive electron and neutrino cross sections. Within our 
Giessen BUU (GiBUU) framework, we aim at a consis- 
tent treatment of the initial vertex and the final state pro- 
cesses and we emphasize that these should not be treated 
separately. 

Compared to our previous work [1^, where only the 
A resonance was included, we have improved our model 
considerably: Besides the A we have included 12 higher 
resonances. Recent electron scattering data were used 
to fix the resonance vector form factors; but also the 
axial ones were refitted. A non-resonant background is 
accounted for also in the neutrino case. Semi-inclusive 
processes such as pion production or nucleon knockout 
will be discussed in a forthcoming publication and are not 
considered in this work. First results, where we included 
only QE and A excitation, have already been presented in 
Refs. [13 for CC and NC neutrino induced reactions. 

This article is organized as follows. First, we present 
our model for the electron and neutrino nucleon interac- 
tion where we put the emphasis on the neutrino-micleon 



reaction. Then, we discuss the nuclear ground state and 
the modifications of the elementary lepton-nucleon ver- 
tex inside the nucleus. Afterwards, results for both elec- 
tron and neutrino inclusive scattering off Oxygen are pre- 
sented. Main attention is given to the electron-nucleus 
reaction where experimental data allow for conclusive 
comparisons. The effect of different nuclear corrections 
is discussed. Finally, we compare our model to other ap- 
proaches. The appendices complete Sec. II by giving all 
necessary details. 



II. LEPTON-NUCLEON INTERACTION 



In this section, we explain in detail our model for 
the lepton-nucleon interaction. In the region of inter- 
mediate lepton beam energies {Eheam ~ 0.5 — 2GeV), 
the cross section contains contributions from quasielas- 
tic (QE) scattering {tN £'N'), resonance (R) excita- 
tion {£N — > £'R) and direct, i.e. non-resonant, single-pion 
production (£N £'ttN') treated in our description as 
background (BG). Thus we assume 



dcTtot = dcTQE + ^ dan + dfiBG ■ 



(1) 



The dynamics of the interaction is encoded in the ab- 
solute value of the matrix element squared, summed and 
averaged over initial and final spins. 



\M 



QE,R,BG| 



EM,CC,NC^M''^QE,R,BG 



(2) 



with C'em = 47ra/g^, Cqc = Gp cosOc/V^ or Cnc = 
Gp/V^. is the four-momentum transfer; a = 1/137, 



Gf ^ 1.16637- 10-5 GeV" 
leptonic tensor is given by 



and cos 9, 



c 



0.9745. The 



1 



-Tr 



+ mi)l^{l(:' + mt')lu 



(3) 



where = 7^(1—075) and = 7o'^7o; ^ (^0 denotes the 
4- vector of the incoming (outgoing) lepton and (m^/ ) 
the corresponding masses. The parameter a depends on 
the reaction process: a = for EM and a = 1 for CC or 
NC neutrino scattering. 

The hadronic currents in i?QE r eg have to be 
parametrized in terms of form factors and thus depend 
not only on the final state but also on the specific process, 
namely electromagnetic (EM) {e~N e~X), charged 
current (CC) (i^gN — > £~ X) or neutral current (NC) 
{vN ^ vX). 
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A. Quasielastic scattering 

The cross section for quasielastic scattering 
i{k)N{p) -> l'{k')N'{p') is given by^ 



do-QE 



Ik'l 



5{p'' 



M' 



I2\ 



duj dVlk' 327r2 



[(fc-p)2-m2M2] 



1/2 



I-^QEP; (4) 



with M — yjp^ and A/' — yp^ . In the case of free 
nudeons, we have M = M' = Mn, where Mjv denotes 
the average nucleon mass for which we take 938 MeV; 

k° - k'" is the energy transfer and Qk' is the sohd 
angle between incoming and outgoing leptons. 

The hadronic tensor H^'^ for quasielastic scattering is 
given by 



1 

-Tr 
2 



with 



^ QE 



QE 



Fq^ has a V 



A Lorentz structure 



^ QE 



^ "^QE 



-^QE ■ 



with the vector part 



_f2 

2M 



N 



and the axial part 



Mn 



(5) 



(6) 



(7) 



(8) 



(9) 



Here, = p'^ — Pfj.- Pi (« — 1,2) stands either for the 
CC form factors , for the NC form factors or the 
EM form factors Ff^ with N ^p,n; Ta for the CC form 
factor Fa and the NC form factor F^ (analogous for 
Tp). All form factors depend on = —q^. 

The vector form factors F^ and Ff^ can be related 
to the electromagnetic Dirac and Pauli form factors Ff^ 
with N ~ p,n as listed in TableHl Those can be rewritten 
in terms of Sachs form factors, for which we take the 
updated BBBA-2007 parametrization [3l[. The relations 
for the axial form factors are also given in Table [H In the 
cross section formula, J-p appears only multiplied by the 
mass of the outgoing lepton, so it can be safely ignored for 
NC interactions. In the case of CC reactions, we assume 
pion pole dominance, and use the partial conservation of 
the axial current (PCAC) to relate Fa and Fp (details 
are given in Appendix I A 2 a[) . 



MQ') = j^^FAiQ') . 



(10) 



^ Throughout this article we work in the lab frame if not noted 
otherwise. 



For the axial form factor, we assume a standard dipole 
form 



Fa{Q-')^9a[1 + ^^ 



(11) 



with the axial vector coupling constant gA ~ —1.267 [32| . 
The so-called axial mass Ma has recently been refitted by 
Kuzmin et al. [1^ using the BBBA-2007 parametrization 
for the vector form factors. We take their best fit value 



Ma = 0.999 ± 0.011 GeV. 



(12) 



The strangeness content of the nucleon is contained 
in F( 2 and F| (see Table |l| which are present in NC 
neutrino scattering. No definite conclusion can be drawn 
from data, thus we set for simplicity 



j^r(o) = o, 

J^2(0)=0, 



and 



fi(q2 



As 



(13) 
(14) 



(15) 



assuming that the strange axial mass is equal to the non- 
strange one. Here As denotes the strangeness contribu- 
tion to the nucleon spin for which we use As — —0.15. 



B. Excitation of baryon resonances 

This section is devoted to the second part in our gen- 
eral decomposition of the cross section given in Eq. ([1]), 
namely the excitation of resonances Aup. As will be 
shown, the electromagnetic form factors are taken from 
the MAID analysis [SJ, |3^[3^. In this analysis, 13 res- 
onances with invariant masses of less than 2 GeV are 
included. This limits then also the number of resonances 
in our model to 13 — they are listed in Table [III 

The cross section for resonance excitation £{k)N{p) 
e'{k')R{p') is given by 



dgR 
dw drik 



327r2 



Aip'') 



[{k-py 



i2M21^/' 



(16) 



with M = \/p2. A{p'^) denotes the vacuum spectral 
function of the particle, which is given by a Breit-Wigner 
distribution 



Aip") 



I p- 



TT (p'2-M2)2+p'2r2(p/) 



(17) 



with the momentum-dependent width F taken from the 
Manley analysis [s^l (see Table |TT] for the values at the 
pole). 
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reaction 


replace JT^ with 

for J = 1/2 




for 1 = 3/2 


replace J^a with 
for 1 = 1/2 for 1 = 3/2 


C p — ^ 6 












e'ri ^ e-X° 


















FY = ^VSFr 




V3F4 


vn^ rx+ 


pV ^ pp _ pn 




pV - _pN 


Fa 


Fa 


1/p vX'^ 


Ff = (i - 2sin2 eJ)F[ - \Fr - \F! 




= -(1 - 2sm^ ew)Fl' 






vn uX° 


Fr = {h~- 2sm2 ew)Fr ~ \Ff - \F! 




= -(l-2sm^ ew)F!'' 




F^ = Fa 



TABLE I: Isospin relations for the form factors. X stands for the nucleon A'^, and N* or A resonances. J-i is the generalized 
vector form factor in Eq. ^ and Eq. (I22p . which has — depending on the process — to be substituted following the prescription 
in the table (analogous for J^a in Eq. (|9} and Eq. (|23|) ). Note that in the case of an isospin 1/2 — > 3/2 transition, the form 
factors are equal for proton and neutron which is indicated by the index N (instead of p or n). This replacement scheme 
is identical for both, spin 1/2 and spin 3/2 resonance excitations. In the case of spin 3/2, replace JT with C 
(cf. Eq. (|3H) and Eq. (|32p). and F with C. A detailed derivation of these relations is given in Appendix I A 2 al for 1 = 1/2 
resonances and in Appendix I A 2 bl for 1 = 3/2 resonances. 



name 


Mr [GeV] 


J 


1 P 


rf)°* [GeV] 


nN branching ratio 


Fa(0) or C^iO) 


P33(1232) 


1.232 


3/2 


3/2 + 


0.118 


1.00 


1.17 


Pii(1440) 


1.462 


1/2 


1/2 + 


0.391 


0.69 


-0.52 


Di3(1520) 


1.524 


3/2 


1/2 - 


0.124 


0.59 


-2.15 


Sii(1535) 


1.534 


1/2 


1/2 - 


0.151 


0.51 


-0.23 


S3i(1620) 


1.672 


1/2 


3/2 


0.154 


0.09 


0.05 


Sii(1650) 


1.659 


1/2 


1/2 - 


0.173 


0.89 


-0.25 


Di5(1675) 


1.676 


5/2 


1/2 - 


0.159 


0.47 


-1.38 


Fi5(1680) 


1.684 


5/2 


1/2 + 


0.139 


0.70 


-0.43 


D33(1700) 


1.762 


3/2 


3/2 


0.599 


0.14 


0.84 


Pi3(1720) 


1.717 


3/2 


1/2 + 


0.383 


0.13 


-0.29 


F35(1905) 


1.881 


5/2 


3/2 + 


0.327 


0.12 


0.15 


P3i(1910) 


1.882 


1/2 


3/2 + 


0.239 


0.23 


0.08 


F37(1950) 


1.945 


7/2 


3/2 + 


0.300 


0.38 


0.24 



TABLE IL Properties of the resonances included in our model. The Breit-Wigner mass A1_r, spin J, isospin 1, parity P, the 
vacuum total decay width Fg"' at the pole, the branching ratio into ttA and the axial coupling (Fa(0) for spin 1/2 states; C^(0) 
for states with spin 3/2 or higher) are listed. The resonance parameters are taken from the analysis of Mauley et al. [371 ]. 



1. Excitation of spin 1/2 resonances 



and for states with negative parity (e.g. Sii(1535)), 



The hadronic tensor for the excitation of a spin 1/2 
resonance is given by 



-"1/2 



1. 



-Tr 



/2V 



with ill = ^^^'^ = VF^ and 

For states with positive parity (e.g. Pii(1440)), 

^ 1/2+ ~ '^1/2 -^1/2 



(18) 

(19) 
(20) 



^ 1/2- 



'^1/2 -^1/2 



where the vector part V^yj is given by 



(2M^)2 



_f2 

2M. 



N 



and the axial part A^i2 by 



(21) 



(22) 



(23) 
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As in the QE case, JF^ (i = 1, 2) stands either for the CC 
form factors F-^, for the NC form factors Ff^ or the EM 
form factors Ff^ with N — p,n; analogous for !Fa and 
Tp. The form factors depend on = —q^. 

FY can be related to the electromagnetic transition 
form factors Fl^ with N = p,n as listed in Table HI thus 
the form factors for neutrino and electron scattering are 
related. The form factors Ff'^ can be derived from he- 
licity amplitudes extracted from electron scattering ex- 
periments. The explicit relations between the form fac- 
tors Ff'^ and the helicity amplitudes ^^/2' ^'^'^ '^1/2 
given in Appendix IB II for both, positive and negative 
parity states. We use these relations to extract the form 
factors from the results of the recent MAID2005 analy- 
sis [33, m, [3^ HI] for the helicity amplitudes and their 
Q^-dependence. 

Experimental information on the N — R axial form 
factors Fa and Fp is very limited. Goldberger-Treiman 
relations have been derived for the axial couplings [39| . 
but there is no information about the Q^-dependence. 
We will follow this approach and apply PCAC and pion 
pole dominance to derive the axial couplings Fa{Q) and 
to relate Fa and Fp. The derivation, performed in Ap- 
pendix [Cjal leads to 

(4^ + 

with + (— ) for positive (negative) parity resonances. The 
Q^-dependence of the axial form factor is neither fully 
constrained by theory nor by experiment, so we assume 
— as in the QE case — a dipole dependence 

FA{Qi') = FA{Q)(\ + -^ ; (25) 

with M\ = 1 GeV as for the nucleon. The coupling Fa(0) 
can be related to the 7riVi?-coupling as detailed in Ap- 
pendix IC 2 al (non-diagonal Goldberger-Treiman relation) 
and we obtain the results summarized in Table |TT1 

As in the nucleon case, strange form factors can con- 
tribute for isospin 1/2 —> 1/2 transitions (cf. Appendix 
lA 2 a[) . However, the present experimental status does 
not allow any conclusions on the strange transition form 
factors, thus, we neglect them in this work and set F/ 
and F| to zero. 

Several studies have addressed the properties of the 
Pii(1440) as well as its electroma gnet ic and weak pro- 
duction [S im, S S S, Ei, Hill. In the works of 
Refs. j39l. l48j| the form factors of the four lowest lying res- 
onances, P33(1232), Pii(1440), Di3(1520) and Sii(1535), 
are extracted using helicity amplitudes for the vector 
form factors and PCAC for the determination of the axial 
couplings. 

We now compare our form factors to the most recent 
works by Lalakulich et al. 48] and Hernandez et al. [13] . 
To perform a meaningful comparison, we have corrected 



this work Hernandez Lalakulich 



Pii(1440) Ff(0) -1.96 -3.55 1.43 

Fi"(0) 2.26 0.26 -1.43 

Fi^(O) -4.22 -3.81 -2.86 

F^O) -0.46 -0.50 -0.60 

F2"(0) 0.41 0.34 0.60 

Fa^CO) -0.87 -0.84 1.20 

Fa(0) -0.52 -0.63 0.51 

Sii(1535) Ff(0) 0.85 - -1.24 

Fi"(0) -0.02 - 1.24 

Fi^(O) 0.87 - -2.49 

F|(0) 0.46 - -0.66 

F2"(0) -0.35 - 0.66 

F2^(0) 0.82 - -1.32 

Fa{0) -0.23 - 0.21 



TABLE III: Comparison of our form factors for the Pii(1440) 
and the Sii(1535) to the ones compiled by Hernandez 
et al. and Lalakulich et al. 'is] taken at = 0. The 
Sii(1535) resonance is not considered in the work of Hernan- 
dez. We have corrected for different normalizations and sign 
conventions for the axial current. 



for different normalizations^ and sign conventions: both 
use a different sign in the definition of the axial current. 
The result after these corrections is shown in Table [Till 
We agree in all signs with Hernandez et al., but with 
little differences in the numerical values. For the vector 
form factors, these differences can be attributed to the 
different data set which is used, since their Eqs. (28) and 
(29) coincide with our Eqs. (|B5)) and ((B6)l ^. We also 
agree in the expressions for the Goldberger-Treiman re- 
lation (their Eq. (21) vs. our Eq. (jC14p — note the differ- 
ent sign convention); the difference in the axial coupling 
comes from a different value for the 7r7Vi?-coupling. 

Comparing to Lalakulich et al., we find major differ- 
ences. Ignoring a global sign difference, a relative sign 
difference between F^ and Fa remains for the Pii(1440) 
resonance. This is mainly caused by the different sign of 
Ff , which can be attributed — as already pointed out 
by Hernandez et al. [47j — to the extra minus sign in the 
Si/2 amplitude which is missing in the work of Lalakulich 
et al.'^. Furthermore, they perform the calculation in the 



^ which, in both cases, introduces factors (^-'^^w) /2JUjv \ ^j^j^ 
fj, = Mjv -I- Mn in front of the vector form factors 

^ While we include the minus sign in the amplitude explicitly, 
they account for it when comparing to the MAID analysis. 

* If we also do not consider this minus sign, we find reasonable 
agreement. 
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lab frame, while it should be done in the cm frame f38','49'| 
(this affects 51/2, see also the remarks in Appendix [B|) . 

2. Excitation of spin 3/2 resonances 

The excitation of a spin 3/2 final state is described 
within a Rarita-Schwinger formalism where the hadronic 
tensor is given by 



1. 



-"3/2 ~ 2 



with the spin 3/2 projector 

Aa/3 = - (/ + M') 



(26) 



3 M'2 



and 



1 Pa7/3 - PfSla 

3 M' 



- 3/2 = 10^3/2 70 



7/3 



(27) 



(28) 



For states with positive parity as the P33(1232), we 
have 



3/2+ 



'^3/2 



•^3/2 



75 : 



(29) 



and for the negative parity ones (e.g. Di3(1535)), 



^ 3/2- 



'^3/2 



•^3/2 ■ 



(30) 



In terms of form factors, the vector part is given by 



that in the case of isovector transitions, i.e., isospin 1/2 
— > 3/2 transitions, proton and neutron form factors are 
identical. 

The substitutions for the axial form factors are sum- 
marized in Table [D Pion pole dominance and PCAC 
allow us on one side to relate and Cg, and on the 



other side to fix the coupling C^{Q 
IC 2 bl we show that 



0). In Appendix 



Ml 



(33) 



for both parity states. The form factor appears in the 
cross section only multiplied by the mass of the outgoing 
lepton, such that its contribution is negligible in NC and 
rather small even in CC reactions (except for i^r)- 

For isospin 1/2 ^ 1/2 transitions, strange form factors 
can contribute, but are neglected in this work due to the 
lack of experimental information. 

P33(1232). The axial coupling C^(0) is obtained 
using an off-diagonal Goldberger-Treiman relation (see 
Appendix IC2ap . which yields the values given in Ta- 
ble HIl Since one can not constrain C^{Q'^), Cf{Q'^) and 
C5^('3^)/C'5^(0) from theory, these form factors have to 
be extracted from experiment. The available informa- 
tion for such reactions comes mainly from two bubble 
chamber experiments: the 12-foot bubble chamber at Ar- 
gonne (ANL) [H^, HTjand the 7-foot bubble chamber at 
Brookhaven (BNL)[53, HI . 

There have been many earlier attempts to fit these 
form factors, always under the reasonable assumption 
that the vector form factors are precisely known from 
electron scattering. Most of them adopt the Adler 
model [541 where 



IVl pf ivi ^ 



3/2 



^{9"'q-p~q>n 



and the axial part by 



(31) 



and C^{Q^) = 0, 



(34) 



thus, one remains with the still unknown Q^-dependence 
of which has to be extracted from data. A 
parametrization widely used in the analysis of the neu- 
trino experiments |i50|, i51|, iS^, [E^l is 



-^3/2 



'-3 



M, 



N 



iff 



M 



N 



TV 



75 



(32) 



As before, the calligraphic C stands either for the CC 
form factors C^'^ , i = 3, . . . ,6, the electromagnetic tran- 
sition form factors with N = p,n 01 the NC form 
factors CY' ^ as detailed in Table [J Note that current 
conservation implies = for EM transitions. 

Using information from electron scattering, the C^'" 
can be parametrized in the same way as in the previ- 
ous section for spin 1/2 resonances by relating them to 
the MAID helicity amplitudes. The explicit relations be- 
tween the form factors Cf and the helicity amplitudes 
^1/2' ^3/2 ^i'/2 S^^^^'^ Appendix [B2l Note 



Mr 



(35) 



with a = -1.21 and b = 2 GeV^ (5j 



mass fitted to the ANL data of 



with an axial 
0.98 GeV 

The cross section da/dQ^ obtained with this form factor 
parametrization is shown by the dashed line Fig. [T] (la- 
beled with (1)) in comparison to the ANL data. It clearly 
overestimates the data even though the axial form factor 
was originally fitted to exactly those. 

Different parametrizations for the axial form factor 
have been used, e.g. in Ref. who use also a mod- 

ified dipole 



Ml, 



, (36) 
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and take Ma = 1.05 GeV. Using this in combination 
with our MAID based vector form factors one obtains 
the dash-dotted hne in Fig. [1] labeled with (2). Even 
reducing the axial mass to Ma — 0.95 GeV it is not pos- 
sible to get a good description of the ANL data (cf. dot- 
ted curve labeled with (3)) without reducing C^(0). Fi- 
nally, the double-dashed curve denotes the calculation 
performed with the form factor set both vector and axial, 
of Lalakulich et al. [IH^, which also uses the parametriza- 
tion of Eq. ^ for C^(0). Also in this case, no satisfac- 
tory agreement with the data can be reached. Note that 
in this work, the improvement on the vector form factors 
was not applied to refit the axial form factors. 

Comparing these scenarios to the ANL data, we con- 
clude, that improving on the vector form factors without 
readjusting the axial ones results in a worse description 
of the data. A refit of the axial form factor is therefore 
necessary: We rely on the Adler model, and use Eq. ([55]) 
as a starting point for our fit of the Q^-dependence of 
C^. We further assume that PCAC holds, which means 
that we do not take C^(0) as a free parameter. As we 
rely on PCAC for all other resonance couplings (where no 
data are available) we prefer to keep it also here. In ad- 
dition, this coupling was extracted from the BNL data in 
Ref. [5a] and found to be consistent with the PCAC pre- 
diction; also recent lattice results support the off-diagonal 
Goldberger-Treiman relation [57|. We find a = —0.25, 
h = 0.04 GeV^ and M^ = 0.95 GeV as best values which 
is used in the following as our standard parameter set. 
Comparing the solid line in Fig. [1] to the data, one finds 
good agreement for our new parameters. The result for 
the total cross section is shown in Fig. [2] where we com- 
pare to the available data. Also there the agreement with 
the ANL data is satisfactory. 

We stress that we have neglected a non-resonant back- 
ground which is small in the isospin 3/2 channel vp —f 
fi'TT^p. Both, Sato et al. [l^ and Hernandez et al. [l^ 
find within their microscopic models for the non-resonant 
pion background a correction of the order of 10% in this 
channel while the discrepancy caused by the difference 
between the old and the new axial form factor sets is of 
the order of 30%. 

Di3(1520), D33(1700) and Pi3(1720). As the Q^- 
dependence of the axial form factor can be extracted from 
data only for the P33(1232), we assume for the other 
resonances a simple dipole behavior as done for the / = 
1/2 resonances, 

C^{Q^)^C^{0)(^1 + ^^ (37) 

with = 1 GeV. The couplings C^(0) are summarized 
in Table [Til and Cf are expected to be minor and 
thus are neglected. 



^ Here, the agreement at = is better than in the spin 1/2 case; 
cf. the discussion of the differences at the end of Section III B II 



1 I 1 1 1 r 




' ' ' ' ' ' 

0.2 0.4 0.6 0.8 1 

[GeV^] 



FIG. 1: (Color online) Differential cross section da/dQ^ av- 
eraged over the ANL flux taken from [HH] as function of Q"^ 
compared to the ANL data taken from Radecky et al. [5ll |. 
The curves were obtained with different sets of form factors 
as detailed in the text. Our best result is denoted by the solid 
line. To compare with data, we applied an invariant mass cut 
atW < 1.4 GeV. 
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FIG. 2: (Color online) Total u^^p — > jj. Tv^p cross section as a 
function of the neutrino energy compared to the pion produc- 
tion data of of ANL (Refe. jSOi (•), [Hll (■) and ^ (A)) and 
BNL (^3l (><))• One can see that the above shown isospin 
3/2 channel is totally dominated by the A resonance (dashed 
vs. solid line). No cut on the invariant mass is applied. 



Finally, in Table IIVI we compare our form factors for 
the Di3(1520) resonance to the recent compilation of 
Lalakulich et al. [48[ . The agreement is better than in the 
spin 1/2 case which we compared in Table [nH Focussing 
on the dominating form factor C3, we still do not agree 
in the sign of the EM form factors for reasons already 
discussed at the end of the previous section. However, in 
this case no relative sign difference remains for the CC 
form factors as it was in the case of the Pii(1440). 
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this work 


Lalakulich 


Di3(1520) Cf(0) 


-2.70 


2.95 


C3 (0) 


0.28 


— 1.13 




—2.98 


—4.08 


ci(o) 


2.62 


—1.05 








cr(o) 


4.21 


1.51 


Cf(0) 


-1.17 


-0.48 


C5"(0) 


1.96 


-0.17 


^5^(0) 


-3.13 


0.31 


^5^(0) 


-2.15 


-2.10 



TABLE IV: Comparison of our form factors for the D 13 (1520) 
to the ones compiled by Lalakulich et al. 48] taken at — 0. 



3. Resonances with spin > 3/2 

Any formalism describing resonances with spin greater 
than 3/2 is highly complicated [6(|. We thus make as 
simplified assumption that all resonances with spin > 
3/2 can be treated with the spin 3/2 formalism. As we 
will show their contributions are anyway negligible in the 
energy region of interest in this work. 



4- Results 

In Fig. [21 we show the integrated cross section for CC 
(panels (a) and (b)) and NC (panels (c) and (d)) induced 
resonance production on the proton (panels (a) and (c)) 
and on the neutron (panels (b) and (d)). The dominant 
contribution comes from the excitation of the A reso- 
nance (solid line). 



C. Single- TT non-resonant background 

In the preceding section, we have described how we ob- 
tain the needed resonance contributions. Since it is well 
known that photo-nuclear data also contain background 
contributions we now specify their properties. It is es- 
sential to realize that the background contribution must 
be known for all kinematic variables if one is interested 
in calculating also differential cross sections. 

In the case of neutrinos, the need of such a contribution 
is justified by the fact that we underestimate the total 
pion production cross section in the isospin 1/2 chan- 
nels. We know from the comparison to ANL and BNL 
data that the resonance excitation alone does not account 
for all the pion strength (cf. the discussion in our earlier 
work [29]). However, the discrepancy is small compared 



to the total pion production cross section which is dom- 
inated by the isospin 3/2 channel. 

Let us consider the single- tt non-resonant background 
cross section dcreG- It includes vector, axial and also 
interference contributions 



dcTRG = d(T 



BG 



dfj 



V/A 
BG 



A .-.-V 1 J ^noii-V 

= dcrgQ -I- dCTp 



'BG 



(38) 
(39) 



The vector part is fully determined by electron scattering 
data. The axial and the interference term collected under 
the label "non-V" are only present in neutrino scattering 
and will be fitted to the available neutrino data. 



Vector part 



We first discuss the vector part of the non-resonant 
pion cross section, dugQ. Our strategy will be to 
evaluate the unpolarized pion production cross section 
i{k)N{p) i'{k')TT{k^)N{p'), da%^, and to subtract 
afterwards the resonance contribution. We understand 
this subtracted contribution, which then includes also the 
resonance-background interference terms, as a single-7r 
background denoted as dcr^Q. We are actually assuming 
that the resonances do not interfere among themselves or 
with the background, so that 

d^^G 



dLodflh/dQi 



dLudfli-'dni. 



dujdnk'dflk 



The first term of the rhs, the cross section for 
e{k)Nip) ^ f (fc')vr(fc,)iV(p'), da%^, is given by 



ddjdflk'dQi 



Ik'llk. 



[{k-p^ - m\M 



21-1/2 



5127r5 

x%'2-Af'2) l^f^^l^dfcO (41) 



with p' ^k+p~k' -k,,, \k'\ 



raj, and lk„ 



ri^ and with f2fe^ denoting the solid angle of 

the pion. The ^- function 5(p'^ — M'^) eliminates the 
dfc^ integration. The hadronic tensor H'^^ entering the 
matrix element in Eq. (|4ip is written in a form similar to 
that of resonance production 



rrM"^ - try 



with 



(42) 



(43) 



The hadronic vertex can be parametrized in the most 
general way [HH as 
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E 

i=l 



(44) 
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FIG. 3: (Color online) Integrated cross section for CC (panels (a) and (b)) and NC (panels (c) and (d)) induced resonance 
production on the proton (panels (a) and (c)) and on the neutron (panels (b) and (d)). The different lines indicate the different 
resonances. Shown are only the results for the four lowest lying resonances — the sum of the remaining 9 resonances gives the 
double-dashed line. In the case of CC scattering on protons (panel (a)) only the excitation of isospin 3/2 resonances is possible. 



with (in the notation of MAID [62 



= 217^ (^P^q . (fc, - I) - P . g (fc, - 

M!^ = -ij'{rk^-q-^k>^) , 
Mi" = -2i7^(7'"g • P - dPi^) - 2MnM^ , 



(45) 



= -ij^^q^ - q^r) 

and P^ = {p+p')^/2. The so-called invariant ampli- 
, ^6 



tudes A'^'^ , . . . , Aq^, depending both on the probe and 



the reaction channel, are functions of three scalars which 
completely determine all the 4-vectors at the vertex. We 
choose W = ^/s, = —q^q^ and the CM scattering 
angle 6 between q and k^r as such a set of independent 
scalars. 

The second part of Eq. (|40p is obtained in the following 
way: the resonances are assumed to decay isotropically 
in their rest-frame, i.e. 



47r 



and consequently the single resonance contributions are 
given by 



dwdili./drifc. 



d4 



I I R^Ntt ""fc. 



duidfli.' 47r r 



R 



(46) 



The vector part of the resonance cross section has been 
introduced in the previous section and the solid-angle 
transformation is given by (63l | 



|kri(|k,rb'°-|p'|AScos(0^)) 



(47) 



where = /(k^rjp'). 

Electron scattering. In the case of electron scatter- 
ing, the possible channels are 

(48) 
(49) 
(50) 
(51) 





>- 


-> 








e' 








-> e' 


-mr° 


e~ 


- 


-> e' 





For each of these channels, a set of invariant ampli- 
tudes A^^'^^ can be fitte 
parametrization [3^. [ssl \64 



tudes A^^'^^ can be fitted to data. We use the MAID 



10 



Charged current scattering. In the case of charged 
current neutrino scattering, there are three pion produc- 
tion channels, namely 



vn 
vn 



£-7171+ 



(52) 
(53) 
(54) 



Applying isospin relations to relate iq the known 



one gets (for details, see Appendix I A 3| 







^ J^PTT ,EM 


(55) 






^Ptt" ,EM 


(56) 






^mr°,EM ^jT^^pTT^ ,EM 


(57) 



so that the CC vector part is fully determined. 

Neutral current scattering. The following channels 
contribute to neutral current pion production 



vp - 
vp - 
vn - 
vn 



vp-K 
vnn'^ 
vnn^ 
vpn' 



(58) 
(59) 
(60) 
(61) 



Isospin relations required to relate the invariant ampli- 
tudes j^^'^'^'-^ to the known ^^'^■^'^ ^j-g more compli- 
cated than in the CC case and can be found, e.g., in 
Section III of Ref. [11]. 



2. Non-vector part 

In the case of neutrino scattering, an a priori unknown 
non- vector background is present. In principle, one can 
write down a similar expression as Eq. (I44p for the ax- 
ial current. This would lead to a large number of un- 
known amplitudes Ai to be fixed with data. However, 
the scarcity of available data makes this option imprac- 
ticable. 

Recently, microscopic models for the elementary re- 
action have been developed [H, [11] including besides 
non-resonant terms only the A as intermediate resonance 
state. However, as pointed out before, the background 
contribution to the total pion production cross section is 
small because of the dominant isospin 3/2 channel. Thus, 
in the present work, we rather use a simple ansatz for the 
non- vector background expecting that the final results in 
nuclei are not sensitive to the background details. 

Charged current scattering. As already discussed 
in Section Tub 2| we neglect a non-resonant background 
in the isospin 3/2 channel vp — > fj,~'K~^p where it is only 
a small correction. Furthermore, we assume that dugQ 
have the same functional form, i.e.. 



and dag^"-^ 



da 



BG 



da 



BG 



J ^non-V 



(1 



da 



BG 



(62) 



0.4 



0.3 



0.2 



0.1 



E 
o 

00 
CO 

o 



0.3 



0.2 



0.1 



v„ n -» n ji p 




full model 
y resonances 
only A 



v„ n ^ H ji n 




Ev [GeV] 

FIG. 4: (Color online) Total CC pion production cross sec- 
tions for the mixed isospin channels as a function of the neu- 
trino energy compared to the pion production data of of ANL 
(Refs. [HI (•) and [H (■)) and BNL ([HI (x)). The solid 
lines denote the our full result including the non-resonant 
background following Eq. ^ with b""" = 3 and b"'^'^ = 1.5. 
Furthermore, we show the results for pion production only 
through the excitation and the subsequent decay of all res- 
onances (dashed lines) or through the A alone (dash-dotted 
lines). No cut on the invariant mass is applied. 



depends on the process, 
° The data sets for the 



where the global factor 6^ 
vn —^ £~m:'^ or vn —^ i~p'K 
two CC scattering channels off neutrons allow then the 
fit of the two parameters. With b^^ — 3 and 5'"^^ = 1.5 
a reasonable agreement with the ANL data is achieved as 
can be seen from Fig. |4] where the solid line denotes our 
full calculation. This figure emphasizes again the need 
for a non-resonant background in the isospin 1/2 chan- 
nel. Pion production only through the A (dash-dotted 
lines) is not sufficient to describe the data. The inclu- 
sion of higher resonances increases the pion production 
cross section (dashed lines), but still, the non-resonant 
background is required (solid lines). 

Our numbers are in agreement with general isospin 
considerations. The dominant contribution to the non- 
resonant background comes from the nucleon-pole term 
[H, [b^, an isospin 1/2 channel. Using Clebsch-Gordon 
coeflficients, we find for the yields 7r°/7r+ = 2, which 
nicely corresponds to 6^^ /6"'^ = 2. 

Neutral current scattering. Pion data for NC scat- 
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FIG. 5: (Color online) Total NC cross sections for pion pro- 
duction through the excitation and the subsequent decay of 
all included resonances as a function of the neutrino energy. 
Panel (a) shows the cross sections on protons, (b) the ones on 
neutrons. In addition, we show the data of ANL ^6(3] which 
have to be compared with the vn —> upn~ channel (dashed 
line in panel (b)). No cut on the invariant mass is applied. 



tering are even more scarce than for the CC case. The 
available data is compared in Fig. [5] to our results for 
the resonance induced pion production. We find a good 
agreement — dashed line in panel (b) compared to data 
— already without non-resonant background. In view of 
this, we abstain from fitting the non- vector part to this 
data and neglect a NC background. 



III. LEPTON-NUCLEUS INTERACTIONS 
WITHIN THE GIBUU MODEL 

Over the last two decades, the Giessen theory group 
has developed a Boltzmann-Uehling-Uhlenbeck (BUU) 
transport model for nuclear reactions, named Giessen 
BUU (GiBUU) [63, [e^. Besides its original field of ap- 
plication, which were heavy-ion coUisions [H, [2O, [zH, 
the model was also successfully applied to hig h-energy 
non-resonant electron-induced reactions (t^ . [ts! ]. photon- 
, pion- and electron-induced processes in the resonance 
re g i on [tsI . [tgI . [ttI . [ts} , and neutrino- induced reactions 
da ISofT^In April 2008, the GiBUU model source code 
has been published under GNU General Public License 
in a restructured modern form for public use [67j . 



The GiBUU transport approach [67|, |78[ gives a micro- 
scopic description of the final state process. Thereby all 
kind of coupled channel rescattering effects (e.g. charge 
exchange processes, resonance production and decays in 
the medium) are included. All particles are propagated 
in hadronic mean fields and full offshell propagation of 
the hadronic resonances has been implemented. 

In the following, we consider the reactions of electrons 
and neutrinos with nuclei in the impulse approximation. 
Thus, we assume that the elementary projectile scatters 
off one single nucleon bound in the nucleus. In the ini- 
tial step the projectile excites a resonance, generates a 
quasielastic event or a non- resonant background process. 
For exclusive or semi-inclusive reactions the resulting 
particle yield of this initial state must then be propa- 
gated through the nucleus. Such final state interactions 
can then be treated within the transport approach. How- 
ever, in this work we consider only inclusive cross sections 
and postpone the treatment of semi-inclusive reactions to 
a forthcoming publication [79j . 



A. Nuclear ground state 

The phase space density of the nuclcons bound in a 
nucleus is treated within a local Thomas-Fermi (LTF) 
approximation, i.e., at each space point the nucleon mo- 
mentum distribution is given by a Fermi sphere 



r'^(r,p) = e(p^^^(r)-|p|). 



(63) 



with radius p^^ = (Stt^p"''')^/'^. The normalization is 
chosen such that the particle density is retrieved by 



P"'^(r)=.9 / r'^(r,p) ^ 



(27r)3 



(64) 



where g = 2 denotes the spin degeneracy factor. In the 
LTF approximation, the Pauli blocking factor is given by 



PpY(r,p) = l-r'^'(r,p). 



1. Density profiles 



(65) 



The density profiles p{r) are implemented according to 
a parametrization collected in Rcf. ^80*1 : the proton den- 
sity is based on the compilation of Ref. _81] from electron 
scattering; the neutron density is provided by Hartree- 
Fock calculations. 



2. Mean field potentials 

Nucleons (and also resonances and mesons) in the 
medium are exposed to hadronic mean-field potentials. 

The nucleon mean-field potential is defined in the nu- 
cleus rest frame and parametrized according to Welke 
et al. [82] as a sum of a Skyrme term depending only 
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Set A [MeV] 


B [MeV] 


C [MeV] 


r 


A [fm-i] 


#1 -29.3 


57.2 


-63.5 


1.76 


2.13 


#2 -287. 


234. 





1.23 





TABLE V: Parameter sets for the Sk yrm e parametrization of 
Eq. (|66|l (for explicit details see Ref. [S^j). Set #1 is momen- 
tum dependent whereas set #2 is momentum independent. 



on density and a momentum-dependent contribution of 
Yukawa-type interaction 



p{v) 



p{v) 



Po \ Po 

2C f dV 5(/"(r,p') + /nr,p')) 



+ 



PO 



(27r)^ 



1 



(66) 



Different possible parameter sets for the latter potential 
have been fixed by Teis et al. [s^l analyzing nucleon- 
nucleus scattering. The ones being used in this work are 
given in Table [V] parameter set #1 gives a momentum 
dependent potential with a nuclear matter compressibil- 
ity C = 290 MeV, while set #2 results in a momentum 
independent one with C = 215 MeV. Photon-nucleus in- 
teractions indicate that the A(P33(1232))-resonance po- 
tential has a depth of about —30 MeV at po Com- 
paring to a momentum independent nucleon potential, 
which is approximately —50 MeV deep, the A potential 
is, therefore, approximated by 



^A(p,r) 



VN(p,r) 



(67) 



We assume for all spin 3/2 resonances the same potential 
as for the A-resonance. For all spin 1/2 and spin > 3/2 
resonances we assume the same potential as for the nu- 
cleon. Furthermore, we found that the impact of a mean 
field potential acting on the pion (for details cf. Ref. [SGj l 
is negligible for the later regarded processes. 

For convenience, one can rewrite the mean field poten- 
tials as scalar potentials Us which are defined by 



Ml 



y(p,r) = v/p2 + (Afo + [/.(p,r))2 (68) 



with the Breit-Wigner vacuum mass Mq. In this way 
one introduces effective masses for the incoming nucleon 
M = Mm ~\- f7s(p, r) and for the outgoing baryons M' — 
Mb + Usb{p' i'^) where p and p' are the corresponding 
momenta. 



3. Momentum densities 

With a given density parametrization, the single- 
particle phase-space density /(r, p) is fully determined. 



The momentum density is then given by 
n"'^(p)-.9 / r'\r,Y>) 



(27r)3 



(69) 



Fig. [5] shows our result for nP(p) in ^^O as a function of 
the absolute momentum p. To point out the differences 
of the LTF approximation and the common global Fermi 
gas ( GFG) approximation the figure shows also the GFG 
results for comparison. In the GFG approximation one 
assumes that the nucleus is just a sphere with radius i?, 
the one particle phase space densities are given by 

/^/Jr,p) = Q{R-r)Q{pY -p) (70) 

and the momentum densities are simple step-functions 

^l^.^i^f-P)- (71) 



In Fig. El it is shown that the realistic density 
parametrization of Ref. ;8C| leads within the LTF ansatz 
to a less peaked momentum distribution and to more 
strength at low momenta compared to the global Fermi 
gas approximation. Additionally, the LTF approxima- 
tion generates a space-momentum correlation for the 
phase space density, which does not exist for the GFG 
approximation: in Fig. [7] one observes that the proba- 
bility N(j)^ p) to find a nucleon with absolute momen- 
tum j> in a nuclear environment of density p has a ridge 
structure. The low-momentum nucleons experience a low 
density environment while the high momentum nucleons 
tend to sit in a high density environment. This feature 
will become important in the later discussion of inclusive 
electron scattering. 



4-. Spectral functions and in-medium widths 

The spectral function of an outgoing particle with four- 
momentum p is given by 



A{p) 



— ImS(p) 



(72) 



(p2 _ Ml - Rel](p))2 + (ImE(p))2 

where S denotes the self energy. 

To deduce the imaginary part of the self energy of a 
particle in the medium, directly related to its width, we 
have to consider the modification of its free decay width 
due to the Pauli blocking of final state nucleons and also 
collisional broadening due to interactions with the sur- 
rounding nucleons. 

To estimate this collisional broadening of a particle 
with momentum p and energy we employ the low- 
density approximation 

rcoii(Po,P,r) = X! / /'(P'''') f^«(Po,P,p') 

X z;rei(po,P,p') d^p'. (73) 
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FIG. 6: (Color online) Momentum density distribution of pro- 
tons in ^®0. The solid curve represents the LTF result with the 
density parametrization of Ref . [s^ . The dashed and dotted 
curves show the results for a global Fermi momentum of 0.2 
GeV and 0.25 GeV according to Eq. (ITlf) . The normalization 
condition is given by Att f d|p||ppn(|p|) — A. 
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FIG. 7: (Color online) Probability density A'^(jp|, p) of finding 
a nucleon with absolute momentum |p| at a position with 
density p in an ^^O nucleus within the LTF scheme. The 
normalization condition is J dpd|p|A'^(|p|, p) = 1. 



Here the variable Wrei denotes the relative velocity of the 
regarded particle and a nucleon with momentum p'; the 
nucleon phase space densities /* have been introduced 
in Eq. ([70]) . The total nucleon-particle scattering cross 
sections ct^ are chosen according to the GiBUU collision 
term [g^]. Altogether, the imaginary part of the self en- 
ergy is in the rest-frame of the particle given by 

lmY.{po, p, r) = - v^{rfrcc,Pauii blocked (po, P, r) 

+7rcoii(po,P,r)} , (74) 

where 7 denotes the boost factor from nucleus rest frame 
to particle rest frame. The vacuum decay widths Ftioc are 
parametrized according to Manley et al. [s^l (for details 
cf. Ref. [Zl). 

Since the self energy is an analytic function of poj we 



can use dispersion relations to deduce the off-shell behav- 
ior of the real parts from the imaginary part. We apply 
a once-subtracted dispersion relation with the on-shell 
(OS) energy, which is defined by 



„os 



p2+M2-f ReE(pOS^p,r), (75) 



as subtraction point. This yields 

OS 

Rel](p, r) = ReS](pO^ p, r) + 



X p 



ImS(po,P,r) 



-I- Re Co 



(76) 



In the same line as Lehr et al. [871, l88[ , we demand that 
the in-medium shift of the on-shell energy is determined 
by the mean-fields 



- Vp2 + (A/o + f/.(p,r))2, 



(77) 



with the scalar potential Us defined in Eq. ([55]) . Conse- 
quently, the non-dispersive contribution to ReS is given 
by 

ReS(p0^p,r) = 2Afo{/.(p, r) + C/,(p, r)^ . (78) 

In the numerical realization, we approximate the disper- 
sion integral Eq. ([75]) for ReE by 



Rel](p,r) =ReI](p[/'',p,r 




Iml](p[„p,r) 



(P^-P?')(p;,-Po) 
Iml](p;),p,r) E2-p'n \ 

(79) 

with the cutoff parameters Ei — 5 GeV and E2 — 7 GeV. 
Smin is determined by the mass of the lightest decay 
product of the particle. We have checked the dependence 
of our results on the cutoffs and found only a marginal 
impact. The whole procedure guarantees analytical self 
energies and, therefore, normalized spectral functions. 

For the spectral function of the initial state nucleon, we 
consider only the real part of the self energy generated by 
the mean- field potential and neglect the imaginary part. 



B. Lepton- nucleon interaction in the medium 

As outlined in the beginning of this section, the wave- 
length of the exchanged boson is considered small as com- 
pared to intra-nucleon distances in the nucleus such that 
the nuclear reaction may be treated in impulse approxi- 
mation. In this picture the lepton interacts with a single 
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nucleon being embedded in the nuclear medium.^ Hereby 
either a nucleon is scattered in a quasielastic event, or a 
resonance or one direct pion (background contribution) 
are produced. For these reactions we account for the 
in- medium-modifications introduced in Section [ill Al 

In the following we discuss how the different contribu- 
tions to d(Ttot of Eq. Il]) get modified when the nucleon 
is bound inside a nucleus, i.e., 



dfTt,, 



1 ^mcd 



(80) 



where dcr^Q^*^ includes the nuclear modifications discussed 
below. 



1. Quasielastic scattering and resonance excitations 

For the reaction cross sections with a single hadronic 
final state, e.g., resonance excitation or quasielastic scat- 
tering, the implementation of medium-modifications is 
straightforward. The vacuum spectral function for the 
resonance in (|16p and the one of the outgoing nucleon in 
Eq. (gl), namely A^^cip) = - ^n)^ get replaced by 
the in-medium spectral functions 



A V 



A 



medium • 



(81) 



The four-momentum p' of the outgoing hadron is directly 
given by energy and momentum conservation. M and M' 
in the corresponding equations turn into effective masses 
M M + Us due to the presence of the mean field 
potential. 

As detailed by Naus et al. [O^l, the most general in- 
medium vertex for the quasielastic 7*A'^ — > N' process 
has in general 12 linearly independent Lorentz structures, 
which yields also 12 different form factors. These form 
factors may, unlike the vacuum ones, depend on all the 
possible independent Lorentz scalars — e.g. q^, and 
P ■ q where P = p+p'. Unfortunately, there is no feasible 
way to extract all these different form factors from the 
available experimental data. Therefore, we model the 
in-medium nucleon-photon vertex in the same spirit as 
de Forest [Olj, which means that the in- medium vertex 
structure is assumed to be the same as the vacuum one. 
Furthermore, we assume that the vertex structure and 
the form factors are not modified by the medium and 
depend only on Q^. 

In the case of quasielastic scattering, the effective 
masses of in- and outgoing nucleons {M and M') are 
in general not equal due to the momentum dependence 
of the mean-field potential. Adding an extra term in the 
QE current by replacing in Eq. ([8]) 



Q 



(82) 



ensures vector current conservation even when the masses 
of the initial and final nucleons differ. This extra term 
vanishes by applying Gordon identities in the limit of 
free nucleons, so that when Al ~ M' , we retrieve the 
standard expression of Eq. ([8]) . 

Pollock et al. [o^] find that this extra term ^g^/Q^, 
which guarantees charge conservation at the vertex, can 
also be generated via a gauge transformation of the re- 
sult obtained in Landau gauge to Feynman gauge. They 
argue that the above term is no more than a gauge relict. 
We have checked the impact of this term on our results 
finding no noticeable modifications. 

Note that in this model the bound state properties of 
the nucleon are present for both for the initial and final 
states. We also note that a struck nucleon does not nec- 
essarily has to be knocked out of the nucleus, but can also 
stay bound in the potential. This is different from the 
model of Benhar et al. [H, [l^ , where the outgoing nu- 
cleon is always assumed to be knocked out and described 
in terms of a single parameter, an average removal en- 
ergy. 

As in the quasielastic case, we also assume that for res- 
onance production, the in-medium vertex structure and 
the form factors are not modified by the medium. 



2. Single-pion production backgrounds 



Let us now consider the single pion production back- 
ground channel. The single pion cross section obtained 
from MAID is considered as a parametrization of the vac- 
uum data and the background is the difference of data 
and resonance contributions. Thus we first construct ac- 
cording to Eq. the total single pion cross section 
using the MAID input and the contribution of the reso- 
nances to single-pion production using vacuum kinemat- 
ics (no modifications besides Fermi motion, in particu- 
lar no potentials, no Pauli blocking). Then we evalu- 
ate the difference of both cross sections using vacuum 
kinematics — this yields the background cross section 
d(TBG/(da;drife'). 

This background cross section is now assumed not to 
be infiuenced by the potentials, which means that the 
in-medium background cross section is given by 



dcTBG, 



medium 



dcTBG, 



vacuum 



dwdilfc' 



(Pvac,g), (83) 



We note that the impulse approximation is expected to work only 

at high momenta (|p| > 300 Me V) where Pvac = 



M2 
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IV. INCLUSIVE SCATTERING OFF NUCLEI 

The inclusive cross section for scattering off a nucleus 
is given by 



da, 



tot 



(2^ 



<d{pF{r) -p) 



=9 I 

/ nucleus 

^—^daZ^'PMr.p). (84) 

Wrcl K p 

where the kinematical prefactor corrects for the nuclear 
flux. The Pauli blocking factor PpB is given in Eq. ([65)1 . 
The relative velocity between lepton and nucleus is Wrci ~ 
1. Within our GiBUU framework, this integral is solved 
numerically using a Monte Carlo method. 

In the following, we first present our results for electron 
and neutrino induced reactions and compare to other 
models. 



A. Inclusive electron scattering 

Electron scattering off nuclei in the regime of energy 
transfers between 0.1 and 1 GeV^ has been addressed by 
several experiments within the last two decades, for a re- 
cent review cf. Ref. [l| . Comparing the measured nuclear 
cross sections to the nucleon cross sections, several mod- 
ifications could be observed. First of all, nuclear Fermi 
motion leads to a smearing of the peak structures such 
in the quasielastic and A regions. Furthermore, one ob- 
serves a quenching of the spectral strength around the 
quasielastic peak, which has also been interpreted as a 
violation of the Coulomb sum-rule [ol, [ol, [o^ . In con- 
trast to the quenching in the peak region one observed 
an enhancement in the so-called dip-region in between 
quasielastic and A peak. The peak position of the A reso- 
nance was found to be both A and dependent [H, [93] : 
a shift towards lower masses for < 0.1 and a sh ift to- 
wards higher masses for higher Q"^ [93, [sl, [51, Il00l | . 



1. Scattering of Oxygen 

Here the results of our model are compared to data 
and the impact of the most prominent model ingredi- 
ents is investigated. In Fig. [5] we show our results for 
the inclusive reaction ^^0(e^, e^) X a.t a beam energy of 
700 MeV and for different nucleon mean field potentials, 
in-medium changes to the width have been neglected. 
The dash-dotted curve denotes the result without po- 
tentials, including only Fermi motion and Pauli block- 
ing. The presence of a momentum-independent potential 
(dashed curve) does not practically change the QE peak 
{u! = 0. — 0.15 GeV). However, the single-pion region 
(w > 0.2 GeV) is modified. This is due to the effect 
that the A — dominating this region — is less strongly 
bound than the nucleon. As a consequence, more en- 
ergy must be transferred to the nucleon to compensate 
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FIG. 8: (Color online) The inclusive electron cross section 
da/iduidfly) on '''O as a function of the energy transfer uj for 
a beam energy of 0.7 GeV and a scattering angle of ^k' ~ 32°. 
The plot shows the results for different nucleon potentials: no 
potential (dash-dotted line), momentum-independent poten- 
tial (dashed line) and momentum-dependent potential (solid 
line). The calculations do not include in- medium changes of 
the widths. 



the binding. When the momentum-dependent mean field 
is included (solid curve), then the faster (on average) fi- 
nal state nucleons experience a shallower potential than 
the initial state ones. Also the resonance potential gets 
shallower for higher momentum. Therefore, even more 
energy must be transferred by the photon leading to a 
broadening of the QE peak towards a higher energy trans- 
fer oj and to a shift of approximately 8 MeV; also the 
single pion spectrum is slightly shifted towards higher 
energies and broadened. A similar resul t has also been 
obtained within the Walecka model jl0l| . There the nu- 
cleon mass acquires an effective mass M*(r) such that 
the nucleon energy is given by Ep = ^p2 + Af*(r) which 

V{r, p) with the 



can be rewritten as Ep = y^p^ 
momentum-dependent potential 



M - 



V{r,p) = Vp2 + M^(r)- ^/p'^+M . 

For small momenta (|p| ^ M*, \p\ ^ M), one obtains 
a simple quadratic dependence of the potential on the 
momentum 



V{r,p)-^-' 



P 



M -M*{r) 
2MM*{v) 



+ M*{r)-M. 



Rosenfelder [10l| suggests that the value of M* can then 
be used to fit the QE peak. We emphasize however, that 
we do not fit our potential to the electron data since it 
has already been fixed by nucleon- nucleus scattering (83| . 

Fig. [5] shows the comparison o f our mod el to the data 
measured by Anghinolfi et al. jlOOl . Il02l |. The over- 
all agreement with data is good, specially when the in- 
medium widths are taken into account. In particular, the 
description of the QE peak is successful except at the 



16 



lowest beam energy (0.7 GeV). At large energy transfers, 
the lack of double pion production strength leads to an 
underestimation of the data 

To analyze the problem at -Ebcam — 0.7 GeV more 
closely, we compare in Fig. [10] our model to m easure- 
ments performed by O'ConncU tt al. [H, Il03l |. where 
the kinematical constraints are similar to the (-Ebcam = 
0.700 GeV, 6'k' = 32.0°)-run performed by Anghinolfi 
et al. To compare experiments one often defines the so- 
called at the QE-peak. It corresponds to the at 
which the center-of-mass energy at the hadronic vertex 
equals the nucleon mass if one assumes a free nucleon tar- 
get at rest. This parameter depends only on the beam 
energy i?bcam and the electron scattering angle 6'k' and 
is given by 



-peak 



= 2M 



N 



^boam(l 



cos 9[f^i ) 



N 



.(1 



cos( 



'k' 



(85) 



Note that it does not include any in-medium input, but 
it gives a simple estimate of the value at which the 
QE maximum is reached. As for the Anghinolfi data, 
we get also for the O'Connell data [ol, Il03l | good corre- 
spondence for the measurement with higher QgE-pcak 
0.190 GeV^ while we overshoot the QE-peak height for 
the lower Q2^-pcak = 0.106 GeV^ 

In Fig. [TTl we show the contribution of the differ- 
ent production mechanisms to the total electron-nucleus 
cross section including all in-medium modifications. The 
quasielastic, single-7r and 2tt contributions to the initial 
scattering process are shown. One observes that when 
going from low to high beam energies, the importance 
of single-TT and 27r production mechanism gradually in- 
creases, whereas at low energies the quasielastic contri- 
bution dominates. Note that this result does not include 
any FSI of the outgoing particles except for the mean- 
field potential and of the outgoing nucleon and the clas- 
sification into different channels is solely based on the 
initial vertex and not on the final-state multiplicities. 

2. Impact of the initial phase space distribution. 

In the following we compare results for different phase 
space distributions of the target nucleons. As outlined in 
SectionHnAl we assume for the nuclear ground state that 
the positions of the nucleons arc distributed according to 
density parametrizations obtained from low-energy elec- 
tron scattering and Hartree-Fock calculations. The mo- 
menta of the nucleons are distributed according to a local 
Thomas-Fermi (LTF) approximation, i.e. it is assumed 
that at each space-point r the nucleon momenta occupy 
a uniform sphere in momentum space with a radius given 
by the Fermi momentum p/ (r ) . In Section FlII A 3l we have 
already compared the simpler global Fermi Gas model to 
LTF finding quite different momentum densities. 

For a nucleon at rest, QE scattering takes place at a 
given energy transfer (for fixed incoming beam energy 
and scattering angle) 
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FIG. 9: (Color online) The inclusive electron cross section 
dcr/(dajdf2j,;) on ^''O as a function of the energy transfer lu at 
five distinct fixed electron energies (0.7, 0.88, 1.08, 1.2 and 1.5 
GeV) and a scattering angle of S^/ = 32° . The dashed lines 
denote our result, where we include all in-medium modifica- 
tions besides collisional broadening. The solid lines denote 
the full calculation, which includes in- medium ch anges of the 
width. The data are taken from Refs. [lOOl . Il02l | and the pa- 
rameter QqE-peak is evaluated according to Eq. (|55)) . 
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FIG. 10: (Color online) The inclusive electron cross section 
da / {duidQk') on ^^O as a function of the energy transfer oj for 
a beam energy of 0.537 GeV and 0.737 GeV and a scattering 
angle of 6^1 = 37.1° in com parison to the data measured 
by O'Connell et al. [9^ . llOSj . The parameter QqE-pcak is 
evaluated according to Eq, 



QE-pcak 

'2M 



^gca.n(l-COsgkO 
bcam(l - COS 61k') + M 



(86) 



The finite target nucleon momenta within a Fermi gas 
lead to a finite range of possible lo centered roughly 
around w = QqE-pcak/^-^i for which a QE-event can 
be realized. The size of this range is determined by the 
Fermi momentum: the larger the Fermi momentum the 
broader lo range is present. 

To compare the LTF and GFG distributions we use the 
very same physics input, in particular the same potentials 
and widths, varying only the initial r and p distributions 
focussing on the quasielastic peak. 

In panel (a) of Fig. [T^l we show the results for differ- 
ent input distributions neglecting all in-medium modifi- 
cations besides Pauli blocking and Fermi motion. The 
dashed-dotted curve represents the result with our stan- 
dard momentum distribution according to the local- 
Thomas-Fermi ansatz. One observes that an initial dis- 
tribution according to the GFG with a constant Fermi 
momentum of 0.2 GeV yields already a slightly lower and 
broader QE-peak compared to the LTF result. Fig. [S] 
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FIG. 11: (Color online) The inclusive electron cross section 
da/idLudny) on ^''O as a function of the energy transfer uj 
at five distinct fixed electron energies (0.7, 0.88, 1.08, 1.2 and 
1.5 GeV) and a scattering angle of O^^i = 32°. The solid lines 
denote our full result, where we include all in-medium mod- 
ifications and in particu lar in-me dium changes of the width. 
The data are taken from [lOOl . Il02l | . The dashed lines show the 
quasielastic contribution, the dotted ones the single- tt and the 
dash-dotted ones the 2n contribution to the initial scattering 
process. This result does not include any FSI of the outgoing 
particles. The parameter QqE-peak is evaluated according to 
Eq. (113. 
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12: (Color online) The inclusive electron cross section 
da / {duidQk') on ^^O as a function of the energy transfer oj for 
a beam energy of 0.7 GeV and a scattering angle of 6k' = 32° . 
The graphs in panel (a) show the results for a calculation 
where we neglected potentials and in-medium width modifica- 
tions. The results for various assumptions concerning the mo- 
mentum distribution of the target nucleons are shown: Fermi 
gas with Fermi momentum p/ = 0.2 GeV (solid line), Fermi 
gas with Fermi momentum Pf = 0.25 GeV (dashed line), mo- 
mentum distribution according to local Thomas-Fermi (LTF) 
approximation (dash-dotted line). The calculations shown in 
panel (b) include mean field potentials, but also no in-medium 
changes of the widths. The GFG result in panel (b) was ob- 
tained with a Fermi momentum of 212 MeV. 



shows the momentum densities of the target nucleons 
for the LTF scheme (solid line) and the Fermi gas with 
Pf — 200 MeV. Note that both distributions have the 
same average for the absolute value of the nucleon mo- 
mentum (150 MeV). However, within the LTF ansatz 
there are more nucleons with low momenta which leads 
to more strength around uj « QgE-pcak/^-^- The surplus 
of high-momentum nucleons for the LTF case compared 
to the GFG leads to a slightly higher maximal photon 
energy for LTF and to a slight broadening of the peak 
at low and high photon energies for LTF. With a further 
increase of the Fermi momentum to 0.25 GeV we obtain 



a further broadening of the QE peak and a prominent re- 
duction of the peak height. The peak position is however 
insensitive to the magnitude of the Fermi momentum, so 
that for a GFG the height of the QE-peak could be fitted 
by a variation the Fermi momentum parameter. 

In panel (b) of Fig. [T^l we additionally consider the 
impact of two different mean-field potentials for the nu- 
cleons and compare the LTF results to GFG results with 
Pf = 212 MeV (this Fermi momentum gives a nuclear 
density of 0.84 fm"'^ which is approximately the same as 
the average nuclear density of an nucleus) . In a GFG, 
the density is constant and, therefore, all nucleons are 
bound by a similar potential which only differs due to its 
momentum-dependence. In the LTF description, however, 
the nucleons feel quite different potentials depending on 
their position since the nuclear density becomes position 
dependent. The solid and dash-dotted curves in panel (b) 
of Fig. [12] show the results for a momentum- independent 
potential. The comparison of our results with such a po- 
tential to those without potential (panel (a)) shows that 
a momentum independent potential has almost no visi- 
ble impact on the results because the potentials for the 
incoming and outgoing nucleons are the same and the 
energy transfer stays constant. Also the slightly lower 
nucleon in-medium masses due to the mean-field poten- 
tial do not lead to a sizable modification of the cross 
section. 

The dashed and dashed-dotted curves on panel (b) of 
Fig. [H] show the results for the GFG and LTF distribu- 
tions when including a momentum dependent potential. 
One observes in both cases a shift of strength towards 
higher photon energies. This shift comes from the po- 
tential difference of the slow and strongly bound target 
nucleon to the faster and less strongly bound final-state 
nucleon. However, there is a qualitative difference in the 
spectra for LTF and GFG: the QE-peak for the GFG 
is shifted by roughly 25 MeV towards higher to while it 
hardly shifts for the LTF ansatz. This feature is caused 
by the fact that the contribution of the slow target nu- 
cleons dominates the region around the in-medium QE- 
peak. Within the LTF scheme, these slow target nucle- 
ons tend to sit at low densities (cf. Fig. [7]) and, there- 
fore, the region around the QE-peak is hardly affected 
by the medium. For the GFG, however, such correlation 
of density and momentum is absent and, both peak and 
off-peak strength are shifted towards higher masses. 

This analysis has shown an important difference in be- 
tween the global Fermi gas and the local Thomas-Fermi 
approximation. We emphasize that in the LTF case, in 
contrast to the GFG, nucleon position and nucleon mo- 
mentum are correlated. If one uses the same momentum- 
dependent potential in both cases, the peak positions of 



Note that a Fermi momentum of 0.25 GeV is not very realistic: 



a realistic density profile for the dilute Oxygen nucleus ISCt l8l| 
gives within LTF only a density-averaged Fermi momentum of 
ca. 200 MeV. 
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the quasielastic peak and also of the resonance peaks 
barely shift for the LTF distribution while there are large 
shifts in the GFG. 



B. Inclusive neutrino scattering 

In Fig. [T3] we show results for CC scattering off 
oxygen. The dash-dotted lines show our result without 
any mean-field potential while it is included in the cal- 
culation denoted with the dashed line. The solid line 
includes, in addition, the in-medium spectral function 
introduced in Eq. (172]) , which accounts for the coUisional 
broadening of the baryons in the medium. The effects 
of these in-medium corrections are identical to the elec- 
tron scattering case shown in Fig. [5] and Fig. [HI We find 
a broadening and a shift of the QE and resonace peaks 
as a consequence of the momentum-dependent potential 
and the in-medium width. Thus, the good description of 
the inclusive electron scattering serves as a direct bench- 
mark for our neutrino calculations even though there are 
still no neutrino data to compare with. 

First results for CC inclusive neutrino scattering off 
nuclei were already presented in Refs. [1^, [2^. There, 
only the excitation of the A resonance was taken into 
account — in this work, we improved on this and include 
13 resonances and a non-resonant pion background. The 
strength of those contributions can be seen in Fig. [141 
The left peaks have their origin in QE scattering (com- 
pare the total yield given by the solid lines to the QE 
yield given by the dashed ones). With increasing lo, one 
enters the pion-production region, which is dominated 
by initial A excitation (solid vs. dash-dotted lines). The 
non-resonant single pion background also contributes in 
this region (dotted lines) while the impact of the higher 
resonances is only visible at high beam energies (double 
dashed lines in lowest panel with Ey = 1.5 GeV). Fig.fTSl 
shows the corresponding plot for NC scattering off Oxy- 
gen. 



C. Comparison to other models 

Finally, we wish to compare our results to the ones ob- 
tained in other models, focusing on the results obtained 
for Oxygen nuclei. The work of Butkevich et al. [l^, 
addresses the quasielastic peak region in a relativistic dis- 
torted wave impulse approximation (RDWIA) approach. 
Their latest model for the ground state spectral func- 
tion includes both a 75% shell-model and a 25% high- 
momentum contribution. When calculating inclusive 
cross sections, the final state wave functions are obtained 
using a real optical potential. This procedure neglects 
the broadening of the outgoing nucleon in the medium 
but incorporates a shift of the outgoing nucleon energy. 
A direct comparison to models based on DWIA with 
complex optical potentials {RDWIA exclusive result in 
Butkevich et al. [2^ or, e.g., Martinez et al. [1^, Maieron 
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FIG. 13: (Color online) Inclusive CC neutrino cross section 
d(T/(dajdf2j./) on ^''O as a function of the energy transfer lu at 
five distinct fixed neutrino energies (0.5, 0.75, 1.0, 1.25 and 
1.5 GeV) and a scattering angle of 0^' ~ 30°. The dash- 
dotted lines denote our calculation without any mean-field 
potential while a momentum-dependent potential is included 
in the calculation shown with the dashed lines. The solid 
lines denote the full calculation, which includes in addition 
in-medium changes of the width. 
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FIG. 14: (Color online) Inclusive CC neutrino cross section 
da / (dujd^lki) on ^''O as a function of the energy transfer cu at 
five distinct fixed neutrino energies (0.5, 0.75, 1.0, 1.25 and 
1.5 GeV) and a scattering angle of = 30°. The solid line 
denotes our full result, where we include all in-medium mod- 
ifications and in particular in-medium changes of the width. 
The dashed hues show the quasielastic contribution, the dot- 
ted ones the single-7r background and the dash-dotted ones 
the contribution coming from an initial A excitation. Higher 
resonance contributions are indicated by the double-dashed 
lines. 
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FIG. 15: (Color online) Inclusive NC neutrino cross section 
da/{dujdQk') on ^^O as a function of the energy transfer lo 
at a distinct fixed neutrino energies of 1 GeV and a scatter- 
ing angle of 6^/ = 30°. The solid line denotes our full result, 
where we include all in-medium modifications and in partic- 
ular in-medium changes of the width. The dashed lines show 
the quasielastic contribution, the dash-dotted ones the con- 
tribution coming from an initial A excitation and the double- 
dashed ones the contributions from higher resonances. 

et al. [l04| ) is not possible since the imaginary part leads 
to a flux reduction in a particular channel. Thus, those 
models rather describe exclusive quasifree single-nucleon 
knockout processes than fully inclusive scattering which 
we consider in this work. 

Meucci et al. [2^, '23*1 apply a relativistic Green's func- 
tion approach to both inclusive and exclusive processes. 
They achieve an impressive description for the data on 
both the longitudinal and transverse response functions 
in the QE region. In the Green's function framework 
a complex optical potential can be incorporated when 
calculating inclusive cross sections without having the 
DWIA-problem of flux reduction. 

The mode l of Benhar, Nakamura and collaborators 
[H, [l^ llOSl . Il06l | is based on non-relativistic nuclear 
many body theory and impulse approximation. It in- 
cludes realistic spectral functions for the hole states ob- 
tained from (e, e'p) data combined with theoretical nu- 
clear matter calculations using the local density approx- 
imation. For the final state spectral function, they rely 
on a correlated Glauber approximation which leads to an 
energy shift of the cross section and to a redistribution 
of the strength (quenching of the peak and enhancement 
of the tail). 

The framework of Benhar et al. provides a state-of-the- 
art description for the hole spectral functions including 
both real and imaginary parts of the self energy. In our 
model, only the real part of the hole self energy is taken 
into account for the initial nucleons (through the mean- 
fleld potential) while the imaginary part is neglected. 

The model of Gil, Nieves and others, applicable to both 
inclusive electron [lj| and neutrino reactions [20], takes 
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into account nuclear corrections beyond Pauli blocking: 
They include long and short range correlations as well 
as particle and hole spectral functions. As in our ap- 
proach, they neglect the imaginary part of the self en- 
ergy in the hole spectral function. Within their model 
Gil et al. 14] achieve very good agreement with data 
for electron-nucleus scattering. Nuclear correlations are 
also taken int o account by Singh, Oset and collabora- 
tors jlOTl . llOSj l for neutrino scattering. They renormalize 
the weak transition strength and are found to be very 
important at low momentum transfers. 

Fig. 1161 shows a comparison of our result including in- 
medium spectral functions (solid line) for inclusive elec- 
tron scattering off ^^O to the results of other models. 
The dotted curves show the latest results of Nakamura, 
Benhar et al. [l^. Their model gives a better descrip- 
tion of the QE-peak than ours at 880 MeV. However, 
at 1200 MeV both achieve equally good descriptions of 
the data in the QE peak. In the pion-production region, 
Nakamura tends to overshoot the data and the A-peak 
position seems somewhat low in w. The data are also 
underestimated in the model of Ref. ^lOSj (dash-dotted 
line). Our model lacks strength in the dip region but 
describes properly the magnitude and position of the A 
peak. The framework of Butkevich et al. [25] , which only 
includes QE scattering, fails in the same kinematical situ- 
ation as our model (700 MeV, double-dashed Hue in upper 
panel) but leads to very good results at 880 MeV. The 
Green's function approach of Meucci et al. [2^ is able to 
describe very well all the data in the QE region. 

In Fig.[T7]we compare our CC calculation to the models 
introduced above for the integrated inclusive cross section 
on Oxygen as a function of the neutrino energy. Panel (a) 
shows the contribution from QE scattering while panel 
(b) shows the pure A contribution. Our full result is 
denoted by the solid line labeled "GiBUU". The over- 
all agreement with the other models is satisfactory. Fo- 
cussing on the QE contribution, our curve is higher than 
other calculations at lower neutrino energies. Also, for 
the A our calculation is slightly higher. However, note 
that differences are already present at the nucleon level 
and do not only arise from the different treatment of nu- 
clear corrections. 

Comparing to electron data, we conclude that our 
model is able to describe the QE region for beam en- 
ergies above 1 GeV while it fails for lower beam energies 
where nuclear many-body effects become important and 
impulse approximation breaks down [89] . From the good 
agreement at higher energies at the QE peak and also in 
the single- TT region, we can conclude, that our low-density 
ansatz for the in- medium width (cf. Eq. (|74p ) and the in- 
clusion of a proper potential incorporate the main nuclear 
corrections. 
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FIG. 16: (Color online) The inclusive electron cross section 
d(j/(dwdnfc/) on as a function of the energy transfer uj 
at four distinct fixed electron energies (0.7, 0.88, 1.08 and 
1.2) and a scattering angle of 0^' = 32°. The solid lines 
denote our full result, the dashed result is taken from Meucci 
et al. [2^, the double dashed one from Butkevich et al. [2^ 
(RDWIA), the dotted one from Nakamu ra e t al. [l^ and the 
dashed-dotted one from Benhar et al. [lOSi ]. The data are 
taken from ^100. . .102.] . 
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FIG. 17: (Color online) Inclusive CC electron neutrino cross 
section a on ^^O as a function of the neutrino energy and 
integrated over the leptonic variables, (a) shows the QE con- 
tribution, (b) the A contribution. Our results are denoted 
by the solid line. We compare to Benhar et al. (their 
"SF-PB" calculation), Singh et al. 107J and Valverde, Nieves 
and others [l09| ] (their full model). 



V. CONCLUSIONS AND OUTLOOK 

In this work, we have presented a model for elementary 
lepton-nucleon scattering and inclusive scattering off nu- 
clei via EM, CC and NC at intermediate lepton energies, 
i.e., beam energies ranging from 0.5 to 2 GeV. In par- 
ticular, the model for the elementary vertex is improved 
considerably compared to our earlier work in Ref. [28j . 

In this energy regime, the scattering process is domi- 
nated by three contributions: quasielastic scattering, res- 
onance excitation and non-resonant pion production. QE 
scattering and resonance excitation (we include 13 reso- 
nances) are both described with a relativistic formalism 
that incorporates recent form factor parametrizations. 
EM, CC and NC form factors are connected via isospin 
relations. For the nucleon vector form factors we apply 
the latest BBBA-2007 analysis accounting for new elec- 
tron scattering data; the resonance vector form factors 
are based on the recent MAID analysis for the helicity 
amplitudes. The axial couplings are obtained applying 
the PCAC theorem. For the A resonance, the axial form 
factor was refitted to the ANL data. For the non-resonant 



pion background, we have used a technique based on in- 
variant amplitudes taken from MAID, allowing us to in- 
corporate background terms not only for EM, but also 
for CC processes, where the additional non-vector parts 
were fitted to the ANL data for total pion production 
cross sections. 

When scattering leptons off bound nucleons, in- 
medium corrections have to be considered. We use a 
local Thomas-Fermi approximation for the phase space 
density of the bound nucleons based on realistic nuclear 
densities. Pauli blocking is naturally included in this 
way. In the initial QE scattering or resonance excitation 
process, we take into account full in-medium kinemat- 
ics, momentum-dependent mean-field potentials and the 
in-medium spectral functions of the outgoing baryons. 
The parameters of the mean-field potential were indepen- 
dently fixed with nucleon-nucleus scattering data. The 
imaginary part of the self energies entering the spectral 
functions are calculated in a consistent way employing 
the low-density approximation. The correct normaliza- 
tion of the spectral functions is ensured by extracting the 
real part of the self energies from dispersion relations. 

Good agreement to the experimental electron scat- 
tering data on Oxygen at beam energies ranging from 
0.7 — 1.5 GeV is achieved both in the QE and in the pion 
production region. The overall agreement is improved by 
taking into account mean fields and in-medium spectral 
functions in addition to a local Fermi gas momentum 
distribution. The momentum-dependence of the mean 
field reshapes the QE peak considerably due to a target- 
momentum dependent energy loss and, in addition, shuf- 
fles strength into the dip region, which is conventionally 
attributed to 2N excitations; also here the description is 
considerably improved due to the coUisional broadening 
of the QE peak. Good description of data is also achieved 
in the single-pion region, while at higher beam energies, 
the data are underestimated at high photon energies u 
due to the fact that non-resonant 27r-production channels 
have not yet been included. 

Still unsolved is the discrepancy at the QE peak at 
very low of < 0.15 wher e our model overesti mates 
both the Anghinolfi ^M, [lo2 and O'Connell [H, [To^ 
data and does not resemble the peak shift. In view of 
this, we studied the influence of the initial momentum 
distribution of the nucleons. We compared our local 
Thomas-Fermi approach to the widely used global Fermi 
gas approximation where the Fermi momentum does not 
depend on the density, and thus, nucleon momentum and 
position are not correlated as in the local Thomas-Fermi 
ansatz. Qualitative and quantitative difference have been 
found, in particular the shift of the QE peak is easily ob- 
tained in the global Fermi gas picture, while the peak 
hardly shifts in the more realistic local Thomas-Fermi 
picture. 

Taking the in-medium modifications on the electron 
scattering results as a benchmark, we have made predic- 
tions for CC and NC inclusive scattering of Oxygen at 
beam energies ranging from 0.5 up to 1.5 GeV. 
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To summarize, we have presented in this article a con- 
sistent model for EM, CC and NC reactions off nucleons 
and inclusive scattering off nuclei in the energy regime of 
interest for current neutrino oscillation experiments. In- 
medium corrections are found to be of considerable im- 
portance and required for a good description of electron 
scattering data, and, therefore, also for neutrino induced 
reactions. The extension of the present model to semi- 
inclusive processes can be performed with the GiBUU 
transport model and results will be given elsewhere. 



for charged currents and 

for neutral currents which includes the isoscalar strange 
axial current Ag . A'^ ^ and A3 are components of the 
same isospin vector. 

2. Isospin relations for QE scattering and 
resonance excitation 



Acknowledgments 



a. Isospin 1/2 1/2 transition 



We thank L. Tiator (Mainz) for his kind support 
and for providing us a compilation of the MAID am- 
plitudes and O. Lalakulich for many discussions. We 
also thank the GiBUU team for instructive discussions. 
LAR acknowledges financial support from the Seneca 
Foundation. This work was supported by the Deutsche 
Forschungsgcmcinschaft . 



APPENDIX A: SYMMETRIES, CURRENTS AND 
ISOSPIN RELATIONS 

1. Hadronic currents 

At quark level, one can directly relate the electromag- 
netic and the weak vector currents assuming only isospin 
symmetry of the strong interactions. If isospin symme- 
try is a good symmetry of the particular hadronic system, 
one expects the obtained relations to be independent of 
the details of the hadronic structure. We then find the 
following structure for the electromagnetic current 



In general, the current is given by 

J'^ = (J^J2^Jr) = v'^-A^ 



with 



(A6) 

(A7) 
(A8) 



where t is the isospin 1/2 transition operator given by 
the Pauli matrices. 

In a transition between isospin 1/2 states, both, 
isoscalar and isovector parts of the current contribute. 
Therefore, we further define for the hypercharge part 



which enters Eq. ([AT]) and Eq. (jAS]) . and 



(A9) 



(AlO) 
(All) 



7^ - It/^ . 
"'em ~ 2 ^ 



"^3 ■ 



(Al) 



is the third component of the isovector (isospin) cur- 
rent and Vy the isoscalar (hypercharge) current. 
The vector part of the charged current is given by 



(A2) 



and the vector part of the neutral current by 



V^c = (1 - 2 sin^ ew)V^ - sin^ OwV^ - ]^V^ , (A3) 

where VP is the hypercharge, Vg is the strange (both 
isoscalar) and 3 is the isovector current. We further 
assume that the matrix elements of the hadronic cur- 
rents ^'ii'i ^3 ^-^^ the same, being related by isospin 
rotation. This is known as "conserved vector current hy- 
pothesis" (CVC). 

Also for the axial part, we assume the hadronic cur- 
rents to have the same structure 



4 A' 
^CC 



A'l 



(A4) 



for the strange component entering Eq. (|A3|) and 
Eq. (jA5[) . respectively. I2 is the 2x2 unit matrix. 

V^, yl^, ... are given in Eq. ^ and © for QE, in 
Eq. ([22]) and ([23]) for spin 1/2 and in Eq. ^ and ^ 
for spin 3/2, but, with the corresponding form factors: 
FY for V^, Fl for V^, ... 

Combining all that, we obtain for the electromagnetic 
transition matrix element 



N* 



"'em 



N* 



* + 



2 



(A12) 



where incorporates the proton form factors F[ . Anal- 
ogously, one finds 



N* 



""em 



(A13) 
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where incorporates the neutron form factors F^. This 
yields the final relations 

= - V/^ , (A14) 
Vi^-^V^ + Vii. (A15) 

The vector part of the charged current transition ma- 
trix element is 



N* 



n) = (N 



(A16) 



which together with Eq. (jA14p relates EM and CC form 
factors. 

Applying the same procedure to NC, one gets 



N* 



TV* 



P 



(1-2 sin^ ew)V^ - sin^ OwV^- - ^ V5 , 

(A17) 

-(1 - 2sin2 9w)V^ - sin^ OwV!^ - ^Vs . 

(A18) 



Hereby, we relate the NC form factors via Eq. (jA14p and 
Eq. (jA15p to the EM proton and neutron form factors — 
in addition, strange form factors coming from V5 have to 
be considered. 

One can proceed in the same way for the axial part 
and finds 



N 



* + 



^CC 



7V*+ 
+ 



+ iA^ 



(A19) 



while 



N* 



4^ 



P 



N* 



4^ 



(A20) 
(A21) 



2 

which connects the NC axial form factors with the CC 
and the strange axial form factors. 

b. Isospin 1/2 —* 3/2 transition 
The full electroweak current can be cast as 



with 



(A22) 

(A23) 
(A24) 



where is the 1/2 3/2 transition operator 



T. 



Tl±iTl 



hi 
2 



±1 



V2 







-'3 ' 



— m 
2 



= I -mix 



hi 
2 



(A25) 
(A26) 



V and are given in Eq. ^ and ^ for spin 1/2 
and in Eq. (pij) and ([5^ for spin 3/2 resonances in terms 
of the corresponding form factors. 

The transition current between 7=1/2 and 7 = 3/2 
has to be purely isovector, so that 



"'em 



^3 



^0 



' 



(A27) 



which yields for the electromagnetic transition matrix el- 
ement 



(A+|j^^|p) = -y|y 
(a°|j^mI") = 



iiio 

2 2 



3 1 
2 2 



7T 7T = -V^ 



1 



(A28) 

= . 
(A29) 



The prefactor of 



is chosen such that the electro- 



magnetic matrix element carries no isospin factor. The 
minus-sign in front is convention — it is co nsiste nt with 
what is commonly found in the literature [llOj . Note 
that we do not include this minus sign in the current, 
but instead we include it in the isospin 3/2 form factors 
(cf. Table HI). 

For the vector part of the weak charged current, we 
obtain 



and therefore 
(A 



VccIp) 



{A+\V^^\n) 



t 

+1 



1 1 

--11 

2 2 



3 3 
2 2 



1 1 
2-2^1 



3 1 
2 2 



\/3V^T_|i, 
(A30) 

= V3V^, 

(A31) 



(A32) 

This agrees with the standard expression found in the 
literature. 

Analogously for the axial part of the charged current, 

A^^ = A^" + lAt^ = VsAf'Tl^ . (A33) 

For neutral current — again, only the isovector part of 
Eq. (|A3|) and Eq. (jA5|) is present — we have 

j;;^ = (1 - 2 sin^ ew)V3^ - ^3 , (A34) 

and thus 



(A+|j^c|f) = (A''|J^cl") = (l-2 



!snr ew)V''-A^. 

(A35) 
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3. Isospin relations for the non-resonant 
background 

In the case of charged current neutrino scattering, we 
have three pion production channels, namely 



vn 
vn 



(A36) 
(A37) 
(A38) 



The vector contribution of the above reactions is 



related to the electromagnetic current J^-^ via [ill 



and for a photon of momentum q moving along the z axis 
^ '|q|,0,0,g") (longitudinal). (B4) 



Our definition of the helicity amplitudes is the same as 
adopted by MAID (cf. Eq. (24) in Ref. [ll; p = q • 3/q° 
as introduced in Eq. (5) of Ref. [i^l). While A1/2 and 
^3/2 are transverse and Lorentz invariant, S1/2 is frame 
dependent. To be consistent with MAID, the resonance 
rest frame has to be used for the calculation [s^ . We use 
the outcome of the MAID2005 analysis^. Note that the 
helicity amplitudes are taken at the Breit-Wigner mass 
Mji, thus we set M' = Af^ in the following. 



1 7'' 

rEM 


n) + 


(piT 






"'em 


P)- 


(piT^ 


"^EM ' 


1. Spin 1/2 



(p-^" l^cch) = {P^° \Jem\p) - ("^° ^emI 



V2{ 



I 7^ I 

rEM I 



(A39) 



For spin 1/2 resonances, only the A1/2 and 51/2 am- 
plitudes are present. 



APPENDIX B: FORM FACTORS AND HELICITY 
AMPLITUDES 



a. Positive parity 



Here we give details on the connection between the 
electromagnetic resonance form factors and helicity am- 
plitudes which we apply to obtain the form factors from 
the helicity amplitudes provided by MAID [si, [H, ill . 

The helicity amplitudes describe the nucleon-resonance 
transition depending on the polarization of the incoming 
photon and the spins of the bary ons; they can be de- 
fined in various ways [4^ . lll2Llll3l |. We use the following 
notation 



^3/2 



'5'i/2 — - 



27ra 



— (R, J, = 1/2 |e+ J^^l N, J, = -1/2) C 
^ {R, J, = 3/2 |e+ J^^l N, J, = 1/2) C , 



(Bl) 



27ra |q| 



X (i?,J, = 1/2 |eOjgj^|iV, J, = 1/2) C: 



In the case of positive parity, the electromagnetic cur- 
rent Jgjy[ is given by the one defined in Eq. (P^ . We 
obtain 




and 



qp.n 
•^1/2 



na {Mm + MrY + Q2 



Y Mjv Ml - 

{Mr - MmY + 



Mr + M, 



AMrM, 



N 



N 



2M, 



N 



Ff'' 



b. Negative parity 



(B5) 



(B6) 



where kR = {M'^ - M^)/2M' and 

/2 ,^2 o2^2 



{M'' 



Ml 



4M' 



(B2) 



M' — \f]7^ is the resonance mass; J^y^ is the electro- 
magnetic transition current. The phase Q is given by 
the relat ive sign between the ttNN and TrA^i? couplings 
[49llll3j |. which we have taken to be -1-1. 
The photon polarization vectors are given by 



In the case of negative parity, is given in Eq. ([22 
with an additional 7^ (cf. Eq. ([2T|l ) and we find 



^1/2 



27rQ {Me 




'N) 



Ml 



Ml 



AMI 



pp.n 



Mr - Mn ^p,„' 



2Mn 



(B7) 



1 



= T^(0, 1, ±i, 0) (transverse) , (B3) 
v2 



* See Section 2 of Ref. [3a | for a history of MAID and the difference 
between different versions. 
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and 



CP," 
^1/2 



Tva {Mm ~ Mji^ +< 



V Mat 
{Mr + Mn)^ 
AMrM, 



Ml 



M^ 



N 



Mr-M^ 



N 



2Mn 



pp.n 



(B8) 



2. Spin 3/2 

a. Positive parity 

In the case of positive parity, the electromagnetic cur- 
rent is given in Eq. pip , an additional 7^ is needed for 
positive parity. This yields 



^1/2 - 



/ 7ra {Mr~MnY + 
3Mn M2 - M 2, 



Af 2, + MnMr + Q2 



M. 



N 



Mr. 



Ml 



Ml 



(B9) 
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^3/2 



M 



N 



Ml 



Ml 



N 



Mr) 



Ml 



J" AT 



Ml 



CI 



Ml 



J" AT 



Ml 



(BIO) 



and 



CP," 

"^1/2 



Tra {Mr- Mn)^ + 



6M 



N 



Ml 



Ml 



^[{Mr - MnY + g2][(Mfl + Mjv)2 + 



Mi 



cr Ml 



Ml 



Ml 



(Bll) 



b. Negative parity 

In the case of negative parity, the electromagnetic cur- 
rent is defined in Eq. ((3T|) and we get 



h/2 



/ na {Mr + MnY + ^ 
SMat' 



Mat 



MAf Affl + 

Mr 



Ml 



Ml 



Ml 



Ml 



121 



(B12) 



/IP," 
"^3/2 



A/at' 



Ml 



^3 {M^-Mr)-^^ 



A^AT 



M^ 



M2 



(B13) 



and 



/2 



Tra (Affl + AfAr)2-h' 



6AfAr 

v4(Mfl 



Ml 



Ml 



M 



N 



Q-'UMr^MnY 



Af2 



Ml 



■Mi 



eg-" Ml + MI + Q 
Ml 2 



(B14) 



APPENDIX C: DERIVATION OF AXIAL plying PCAC according to which the divergence of the 

COUPLINGS 



In this section, we derive the relations between the 
axial form factors and determine the axial couplings ap- 
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axial current is proportional to the pion mass squared. 
Thus, in the chiral limit the axial current is conserved if 
we assume pion pole dominance, i.e., we assume that the 
pseudoscalar part is dominated by the one-pion exchange 
process where a pion is created at the nucleon-resonance 
vertex and then couples to the lepton pair. Its current is 
given by 



— 



(CI) 



where the pion decay constant is given by /tt — 93 MeV. 
For the hadronic vertex, we need the interaction La- 
grangians which are given in the following. 



1. Resonance interaction Lagrangians and widths 

Here we collect the relativistic Lagrangians used for 
the descriptions of the coupling of baryon resonances to 
nucleons and pions. 



a. Resonances with spin 1/2 

The Ri/2Ntt coupling is described by the pseudovector 
Lagrangian 



d^4>-t'f, (C2) 



where the upper (lower) operator holds for positive (neg- 
ative) parity resonances and t = r (t = T''') for / = 1/2 
(/ — 3/2) resonances. 

From this Lagrangian we deduce the following vertex 
factor 



/ j 7''75 



(C3) 



"1^ I J 

The isospin factor is Ciso = \/2 for 1/2 and Ciso = 
— for / = 3/2 resonances, respectively. 

The coupling / can be obtained from the R1/2 ~> 
partial decay width according to 

' An \m„J Mr 

.(C4) 

where Iqcml is the momentum of the outgoing pion in the 
resonance rest frame 



iQcml - • {^5) 

Ir = 1 for isospin 3/2 and Ir = 3 for isospin 1/2 reso- 
nances, and En is the energy of the outgoing nucleon in 
the resonance rest frame given by 

E. = . (C6) 



The couplings f jm-^ are determined via Eq. (|C4p using 
the resonance properties given in Table [TTl To be consis- 
tent with t he ch oice of C in Eq. (|Bip . f /m-^ has to be 
positive Hi, [Til]. 



b. Resonances with spin 3/2 
The i?3/2-/V7r coupling is described by the Lagrangian 

Cr,,,n. = 4r^^. \^^\ d^ct>-t^ , (C7) 



^TT 75 



where Vf^ is a Rarita-Schwinger J'" = 3/2+ field. The 
upper (lower) operator holds for positive (negative) par- 
ity resonances and and t = r (t = T''') for / = 1/2 
(/ = 3/2) resonances. 

From this Lagrangian we derive for the vertex factor 



/ 14 



(C8) 



"^TT [ 75 J 

with the isospin factor Ciso = V2 for /= 1/2 and Ciso = 
^\J\ for / = 3/2 resonances. 

The decay width needed to extract the coupling / is 
— for resonances with = | — given by 



127r \m-^ I Mr 



(C9) 



Ir ( f \- Em±Mm 3 

i R3,2^ttN - 77^3 [ — ] — iqcml ) 

with iQcinl; Ejsi and Ir as above. 

Using the resonance properties given in Table |lTl we 
obtain the couplings f /mT^. 



2. PCAC relations 

a. Resonances with spin 1/2 

We begin with deriving the relation between Fa and 
Fp . Starting with Eq. (multiplied with 7^ for nega- 
tive parity), we obtain 



FA{MR±MN) + ^q^ 



\ 









(CIO) 

where the upper (lower) operator holds for positive (neg- 
ative) parity. PCAC requires that 







(Cll) 



in the chiral limit and, with the extrapolation to non-zero 
pion masses, wc find 



2Mr 



2^ (Mr ± Mn) Mm ^ ,^2^ 



(C12) 



28 



where always the upper sign is for parity + and the lower 
one for parity — . 

To determine the coupling Fa{0), we calculate the pion 
pole contribution to the current and compare that to the 
axial transition current. Starting from Eq. (jClD and us- 
ing Eq. (IC3|) we find 



/ 



= -u{p') Q30 ^ {Mr±Mn) { \ u{p)^^ 



(C13) 

By comparing S^A^ with the divergence of the axial 
current, given in Eq. (jClOp . we obtain in the limit of 



/ 



(C14) 



With this procedure we calculate the axial couplings 
Fa{0) as summarized in Table HH 



where the upper (lower) operator holds for positive (neg- 
ative) parity. Assuming PCAC we can relate the form 
factors (same for both parity states) 



(C16) 



Analogous to the spin 1/2 case, we derive C^(0) by 
starting from Eq. (jCip and using Eq. (|C8|) 



= -^Mp) Ciso l^' } q'' u{p) 2 2 



V2A 



q^ 



(C17) 

Comparing 9^^^ to 9pA!^y2± (cf- Eq- (10151) 1 yields in 
the hmit — + 



(C18) 



6. Resonances with spin 3/2 



We first derive the relation between and Cq . Start- 
ing with Eq. ((32|) (multiplied with 7^ for positive parity), 
we obtain 



9^^3/2± = iV'a(p')9" 



I4 

75 



u{p) , 



(C15) 



The final results for the axial couplings C^(0) are given 
in Table n] 



c. Resonances with spin > 3/2 

We make again the assumption that resonances with 
spin > 3/2 are described with the corresponding spin 
3/2 formula given in the previous sections. The obtained 
couplings are summarized in Table HIl 
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